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1 .  Executive  Summary 

Two  research  topics  have  been  accomplished  through  the  grant.  During  the  period  of 
March  1995  -  March  1998,  research  was  devoted  to  the  microcracking  of  ceramic  matrix 
composites  (CMCs)  at  elevated  temperature.  Analysis  has  been  focused  on  the  investigation 
of  crack  tip  fields  of  matrix  cracks  between  dissimilar  elastic  and  creeping  materials.  A  brief 
derivation  will  be  shown  in  the  next  section.  Depending  on  the  crack  tip  stress  field  and 
strengths  of  the  fiber  and  the  interface,  the  matrix  crack  shown  in  Fig.  a  may  extend  into  the 
fiber  (fiber  failure)  or  kink  along  the  interface  (interface  failure). 
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Fig.  a  Matrix  crack  in  CMCs.  Fig.  b  Failure  mechanism  diagram  for  CMCs 

at  high  temperature 

Fiber  failure  will  occur  if 

o'f  <7e(d  =  0)  ^  °'f  oe{d  =  0) 

— v  ^  and  *  < 

(7*  ae(e  =  90)  X,  Tr6  (6  =  90) 

Otherwise,  interface  failure  will  prevail.  The  interface  failure  can  be  further  divided  into 
debonding  and  sliding  behavior  according  to 

Debonding 
Sliding 

where  a*  is  the  interface  cohesive  strength; 
cr*  is  the  fiber  strength; 


o]_  <79(e  =  7T/2)  ^ae(fl  =  7r/2) 

cr *  <  od(6  —  0)  an  t*  <  Tr9{d  =  nl2) 

r[_  Trf  (0  =  ^/2)  and  r[_  <  Tr9(e  =  7T/2) 

cr’  a e(0  =  0)  an  <Jg(6  =  n / 2) 


r*  is  the  interface  shear  strength. 

If  the  leading-order  solution  is  used  to  evaluate  the  crack-tip  field,  then  failure 
mechanisms  and  the  values  of  ae  {nil)  I  ae  (0)  and  xr8  ( nil )/  oe  (0)  for  different  creep 
exponent  of  n  are  shown  in  Fig.  b.  The  impact  resulting  from  the  research  may  be  summaried 
as  below: 

(1)  All  of  the  composite’s  mechanical  characteristics  that  govern  structural 
performance  and  life  depend  upon  the  constituent  properties  (fibers,  matrix,  interfaces). 
Because  the  constituents  are  variables,  optimization  of  the  property  profiles  needed  for 
design  and  lifing  become  prohibitively  expensive  if  traditional  empirical  procedures  are  used, 
especially  for  CMCs  under  high  temperature  where  creep  deformation  is  prevalent. 

(2)  The  toughness  of  CMCs  at  high  temperature  can  be  enhanced  by  achieving  crack 
deflection  at  the  interfaces  between  fibers  and  matrix,  rather  than  propagation  into  the  fiber. 

In  this  research,  among  the  constituent  properties,  the  critical  parameters,  creep  exponent  of 
the  matrix  (n)  and  strength  ratios,  are  identified  based  on  mechanics-based  models.  These 
parameters  provide  the  designer  with  the  ability  to  vary  interface  properties  as  desired  to 
attain  such  a  goal.  A  failure  mechanism  diagram  is  established  to  demonstrate  the  propensity 
of  these  failure  mechanisms. 

During  the  period  March  -  September,  1998,  the  work  started  with  research  in  the  area  of 
microfailure  of  composite  laminates  with  the  aim  of  establishing  failure  criteria  in 
composites.  The  study  has  been  focused  on  the  investigation  of  microstress  in  a  model 
composite  laminate.  The  model  composite  laminate  consists  of  well-controlled 
microgeometry  of  boron/epoxy  or  SiC/epoxy  layers  bonded  on  an  epoxy  layer.  The 
composite  is  reinforced  by  hundreds  of  fibers  per  square  inch  of  cross-sectional  area.  Since 
the  steep  stress  gradients  occur  in  the  region  of  few  fibers,  only  eight  fibers  have  been 
analyzed.  The  test  of  the  composite  laminate  under  increasing  ey  is  being  performed  at 
Materials  Lab.  at  WPAFB.  By  varying  the  bond  strength  between  the  fibers  and  matrix, 
different  initial  failure  modes,  such  as  interfacial  debonding,  matrix  cracking,  and  fiber 
breakage  can  be  initiated  under  the  macroscopic  loading.  By  correlating  the  analyses  with  the 
experimental  observations  under  critical  loads,  the  failure  initiation  criteria  of  each  failure 
mode  by  incorporating  the  proposed  failure  hypotheses  by  are  currently  evaluated  and  will  be 
quantitatively  established  based  on  stress-based  criteria. 
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2.  Accomplishments/New  Findings 

(A)  Crack-Tip  Fields  for  Matrix  Cracks  between  Elatsic  fiber  and  Creeping  Matrix 


Abstract 

Ceramic  matrix  composites  (CMC)  are  potentially  designed  for  use  in  the  high  temperature 
environment  because  of  their  high  strength  and  noncatastrophic  failure  characteristics.  For 
ceramics  at  elevated  temperature,  the  materials  exhibit  time-dependent  deformation.  The 
phenomenon  is  further  accelerated  by  an  increase  in  the  stress  or  the  temperature. 
Reinforcement  by  incorporating  high-strength  fibers  is  one  of  the  several  approaches  to 
significantly  improving  the  creep  fracture  resistance  of  the  CMC.  Experimental  studies  in 
mechanical  behavior  of  CMC  [1-4]  have  shown  that  failure  is  preceded  by  matrix  cracking  at 
high  temperatures.  Further,  the  creeping  matrix  causes  stress  transfer  from  the  matrix  to  the 
fiber.  However,  there  is  a  lack  of  any  investigation  in  understanding  of  the  relation  between 
the  matrix  cracking  at  the  micro  level  and  the  overall  mechanical  behavior  of  CMC  at 
elevated  temperature. 


In  this  research  effort,  the  crack  tip  fields  for  a  matrix  crack  are  considered.  Using 
generalized  expansions  at  the  crack  tip  in  each  region  and  matching  the  stresses  and 
displacements  across  the  interface  in  an  asymptotic  sense,  a  series  asymptotic  solution  is 
constructed  for  the  stresses  and  strain  rates  near  the  crack  tip.  It  is  found  that  the  stress 
singularities,  to  the  leading  order,  are  the  same  in  each  material,  oscillatory  higher-order 
terms  exist  in  both  regions,  and  stress  higher-order  term  with  the  order  of  0( r°)  appears  in  the 
elastic  material.  The  stress  exponents  and  the  angular  distributions  for  singular  terms  and 
higher  order  terms  are  obtained  for  different  creep  exponents.  A  full  agreement  between 
asymptotic  solutions  and  the  full-field  finite  element  solutions  has  been  obtained. 


Asymptotic  Analysis 

Crack  tip  fields  of  a  matrix  crack  in  an  elastic  power-law  creep  material  perpendicular  to  an 
interface  separating  elastic  material  shown  in  Fig.  1  are  investigated.  The  constitutive 
equation  for  elastic  material  1  can  be  expressed  by 


1  +  v. 


£u 


-Sij+- 


1-2  v, 
3E, 


1 


The  constitutive  equation  for  the  elastic  power-law  creeping  material  2  is  given  by 


e<j  = 


1  +  v. 


1_2v2  3  . 

~S  H - (7,,o  H — £r 

,J  3£^  kk  ,J  2  C 


^0  ) 


/l-l 
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(2) 
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where  e0  and  n  >  1  are  material  constants;  o0  is  a  reference  stress.  Note  that£0  is  allowed  to  be 
a  function  of  time  t  in  modeling  transient  creep  or  material  aging. 
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Fig.  1  Geometry  of  the  cracked  unit  cell 

Referring  to  a  polar  coordinate  system  with  origin  centered  at  the  crack  tip,  as  r  — >  0,  in 
creeping  material  we  seek  an  asymptotic  expansion  of  the  stress  function  02  in  a  form 

<Mg0;O  =  Ao(O^+Vo(^9o)  +  ERe[4(Or‘7‘+2/,(0;^,)]  (3) 

*=1 

where 

<?o  <  Re(<7i )  <  Re(^2 )  <  •  •  • 

To  ensure  the  existence  of  the  solution  in  eq.  (3),  qg  must  be  real  and  qk  (k  >  1)  may  be 
complex.  Thus/,  and  A„  are  real,/  and  Ak  may  be  complex.  For  plane  strain  problem,  the 
stress  component  can  not  be  expressed  directly  in  terms  of  0 .  we  further  assume  the 
deviatoric  stress  szz( r,  6;t)  can  be  expanded  as 

=  V*40)  W  +  X  Re[4  rqt  (0)]  (4) 

i=i 

The  corresponding  expansions  of  in-plane  stress  tensor  a ,  deviatoric  stress  tensors,  stress 
rate  a ,  the  effective  stress  cr  ,  strain  rate  e ,  and  displacement  rate  u  can  be  obtained. 
Therefore,  the  expansions  of  the  compatibility  equation,  plane  strain  condition,  £__  =  0 ,  can 

be  expressed  in  terms  of fk  and  s(J} . 


In  the  region  of  linear  elastic  material,  the  expansion  of  stress  function  0i  is  written  as 


O/r,0;O  =  0o(O£,o+2Fo(0;to)  +  XRe[fi,(O£'‘+2^.(0;/.)]  (5) 

*=i 

The  solutions  for  symmetric  loading  by  satisfying  a  compatibility  equation  are 

Fk  =  ck  cos tk6  +  dk  cos (tk  +2)6  (6) 

where  ck  and  dk  are  constants  to  be  determined. 

The  boundary  conditions  for  this  problem  are 

(°ee)2=(°le])2  =0>  at  6  =  n  (7) 

The  continuity  conditions  along  the  interface,  6  =  Till,  are 

(Oj-K*),  =0.  Ke)2-Ke),  =0*  (»)2-(«),  =0  (8) 
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Asymptotic  Fields 

Collecting  the  coefficients  of  like  power  of  r  in  the  expansions  of  compatibility  equation,  e„ 
=  0 ,  boundary  conditions  and  continuity  conditions,  equating  the  coefficient  of  each  power  of 
r  to  zero,  the  governing  equations  for  each  order  can  be  obtained.  It  is  found  that  the  leading 
order  term  has  identical  stress  exponent,  q0  =  t0,  in  both  regions,  and  higher-order  solutions 
can  be  categorized  into  two  distinct  types.  The  first  type  solution  is  determined  by  a 
homogeneous  ordinary  differential  equation  with  homogeneous  boundary  condition  on  the 
interface,  which  leads  to  linear  eigenvalue  problems.  The  other  is  governed  by  the 
nonhomogeneous  equations.  Comparing  and  then  selecting  the  minimum  value  of  the  stress 
exponents  decided  by  the  eigenvalue  (homogeneuos)  problem  and  by  the  nonhomogeneous 
problems,  we  can  determine  qn  t ,  q„  t:, ...,  etc.  Fig.  2(a)  shows  variations  of  qi  (i  =  0,  1,  2,  3) 
with  n  for  1  <  n  <  20.  By  inserting  the  zero  stress  exponent  in  the  figure,  the  variations  of 
(k  =  0,  1,2,  3, 4)  can  be  directly  obtained. 

Solving  the  governing  equations,  the  stress  exponents  and  angular  distribution  for  the  leading 
order  and  higher  order  terms  can  be  obtained.  The  multi-term  expansion  of  the  crack-tip 
stress  fields  are  given.  For  example,  for  1 .5  <  n  <  2 

(cr)2  =  A0r%ow  +  Axr{2~n)q°  crm  +  r°'2n)q°{A~' Af<y(2)I  +  cT1  Ag-'1  Axoi2)" ) 

+  a~'A'0~n  B2r('~n)q°  0O) 

(a),  =  V9o<t(0)  +  A/2-',)9o(7,1)  +  B2o{2)  +r(3-2',)<7»(A0-,A12(T(3)/  -Hr'A^'V01") 

+  a~'  A'0~"  B2r°~")q°(T{4} 

where  A<Jt)  and  B2(t)  are  real  arbitrary  amplitudes.  A,  =  Al0~"  A0  /  a . 


Discussion  and  Conclusions 

A  unit  cell  in  fiber  reinforced  composite  with  matrix  cracks  under  plane  strain  shown  in  Fig.  1 
has  been  analyzed.  Uniform  displacement  is  applied  to  the  end  of  the  cell  with  geometry  of 
d/a  =  4  and  l/d  =  5  parallel  to  the  fiber  axis  at  t  >  0.  Under  symmetric  loading,  the  boundary 
conditions  for  the  model  are: 


«v  (t)  =  ±le  ( t ) ,  =0  at  y  =  ±l 

ux=0,  0^=0  at  x  =  -a,d/ 2 


where 

£(0  =  <Ul-[l  + 


(n  1  )cc  ^ 


f]  ‘  e„  = 


EfVf 


a  = 


V 


V. 


( 


\  Ym  ) 


EfVf 


\ 


\ 1  +  Em  (\-Vf)j 


Vf  is  the  volume  fraction  of  the  composite.  The  matrix  material  parameters  are  E  =  100  GPa, 
v  =  0.2,  n  =  1.8,  oJE  =  0.002,  eo  =  0.002  /  hr  and  Ob  =  <rc .  The  fiber  material  constants  are  E 

=  350  GPa,  v  =  0.2.  In  this  case,  the  stress  exponents  are  given  by 

q0=~ 0-2940,  ql=(2-n)q0,  q2=(3-2n)q0,  q2=(l-n)q0 
tfi  q^  (2  n)q0 ,  t2  —  0,  (3  2n^q^ ,  t ^  (1  n)qQ 

Numerical  computation  is  implemented  in  the  ABAQUS  finite  element  code.  Eight-node 
isoparametric  elements  with  2x2  Gauss  integration  points  are  used.  The  size  of  the  smallest 
crack  tip  elements  is  1.6  x  10'5a.  Since  the  time-dependent  amplitudes  in  the  asymptotic 
expansion  are  unknown,  the  values  of  the  amplitudes  are  determined  by  matching  asymptotic 
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solution  with  finite  element  solution.  Once  the  values  are  determined,  the  entire  asymptotic 
crack-tip  fields  are  known.  The  matching  is  forced  or  or/and  o$  at  0  =  0  or  7t/2  in  the  fiber. 
Results  show  that  the  values  of  the  unknown  amplitudes  determined  from  different  stress 
components  at  various  locations  vary  slightly. 

Fig.  2(b)  indicates  the  radial  distribution  of  Gee  along  0  =  0  at  t  =  0.06/  £(> .  Here  <7e  is 

normalized  by  Oo-  Figs.  2(c)  and  2(d)  show  that  the  angular  variations  at  r/a  =  2.312  x  10'3 
and  t  =  0.06/  eo .  When  t  =  0.2/  e„  and  r/a  =  1.8  x  10‘3,  the  angular  distributions  are  shown  in 
Figs.  2(e)  and  2(f).  Good  agreement  is  presented  in  these  figures.  Comparison  of  the 
asymptotic  solutions  to  finite  element  solution  indicates  that  the  asymptotic  solutions  can 
describe  the  crack-tip  field  over  a  significant  scale. 

Analysis  demonstrates  that  for  short  time  after  application  of  the  loading,  several  terms  are 
needed  to  characterize  the  crack-tip  fields.  As  t  — >  °°,  stress  field  in  elastic-creeping  material 
can  be  characterized  by  the  leading-order  singular  term  within  a  substantial  distance  from  the 
crack- tip;  while  in  elastic  material  two  parameters  including  the  amplitude  of  the  constant 
stress  term  may  be  needed. 

In  this  paper,  the  generalized  separable  solutions  for  the  stress  functions  <£>,  are  used  which 
lead  to  multi-term  expansion  of  crack-tip  fields.  Although  the  leading  order  solutions  do  not 
have  oscillatory  nature,  the  oscillatory  behavior  may  prevail  in  the  higher-order  solutions. 
Asymptotic  analyses  indicate  that  the  leading  term  is  characterized  by  a  real  parameter  A0(t). 
In  the  creeping  material,  when  n  <7,  the  second-order  term  is  controlled  by  both  elasticity 
and  creep  and  characterized  by  A0\  when  n  >  7,  the  complex  stress  exponent  appears  in  the 
second  term  which  is  described  by  a  complex  parameter.  In  the  elastic  material,  there  exists 
constant  stress  field  with  order  r°  ,  this  term  constitutes  the  second  order  solution  for  the 
range  2  <  n;  the  complex  stress  exponent  enters  the  third  term  for  7  <  n  <  20.  The  analysis 
can  be  extended  to  the  cracks  on  the  interfaces  between  two  elastic-power-law  creep 
materials  (n,  >  n2  >  1). 

References 

1.  R.  T.  Bhatt  and  R.  E.  Philips,  Laminate  Behavior  of  SiC  Fiber  Reinforced  Reaction 
Bonded  Silicon  Nitride  Composite,  NASA  TM  101350,  1988. 

2.  J.  W.  Holmes,  Tensile  Creep  Behavior  of  a  Hot-Presses  SiC  Fiber-Reinforced  SLN4 
Composite.  J.  Materials  Science,  1991,  26:  1801-1814. 

3.  D.  A.  Jablonski  and  R.  T.  Bhatt,  High-Temperature  Tensile  properties  of  Fiber 
Reinforced  Reaction  Bonded  Silicon  Nitride.  J.  Composites  Technology  &  Research, 
1990,  12:  139-146. 

4.  T.  Mah  et  al.,  High-Temperature  Mechanical  Behavior  of  Fiber-Reinforced  Glass 
Ceramic  Matrix  Composites,  J.  American  Ceramic  Society,  1985,  68:  C248-251. 


6 


Fig.  2  (a)  Variation  of  the  stress  exponents  qj  with  n  in  elastic-creeping  material 

(b)  Radial  distributions  of  stress  OQQat  0=0  and  t=0.06/  e0 

(c)  and  (d)  Angular  distributions  of  stress  Gqq  at  r/a=2.312  x  10'3  and  t=0.06/eo 
(e)  and  (f)  Angular  distributions  of  stress  Gqq  at  r/a=1.8  x  10'3  and  t  =  0.2/  e0 
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(B)  Re-examination  of  Microstress  in  a  Model  Composite  Laminate 

Abstract 

Failure  modes  of  composite  laminates  have  their  origins  in  the 
micromechanical  domain.  The  analysis  of  failure  in  fiber  composites  has  traditionally 
followed  two  distinct  levels  of  abstraction.  The  areas  of  investigation  are  known  as 
micromechanics  and  macromechanics: 

The  micromechanics  approach  aims  at  the  involvement  of  microscopic 
inhomogeneities  in  various  kinds  of  micro-failure  processes  by  taking  the  composite 
microstructure  into  account.  The  advantage  of  the  micromechanics  representation  is  that 
detailed  information  is  directly  obtained  about  the  local  interaction  between  the  constituents 
and  micro-failure  mechanisms  (Yuan  and  Yang,  1997).  The  numerical  modeling  of 
exceedingly  complicated  geometric  detail  of  each  fibers  and  matrix,  however,  often  requires 
exceedingly  fine  grids  and  hence  results  in  excessive  computer  cost.  It  is  obvious  that 
conducting  a  stress  analysis  for  composite  laminates  using  this  approach  is  an  almost 
impossible  task  even  beyond  the  computational  capacity  of  the  latest  supercomputers.  Hence, 
the  micromechanical  model  is  mainly  restricted  to  the  strength  prediction  in  the  lamina  level 
or  unidirectional  composites. 

Since  the  inclusion  of  all  inhomogeneties  in  a  mechanical  model  of  laminates  is  an 
overwhelming  task,  in  the  macromechanics  approach,  the  details  of  the  composite 
microstructure  are  ignored  and  the  composite  is  treated  as  a  “locally”  homogeneous  (it  may 
be  inhomogeneous  in  the  laminate  level)  anisotropic  material.  The  homogeneous  anisotropic 
constitutive  relationships  are  determined  from  micromechanics  analysis  of  homogenization 
via  continuum  mechanics.  In  this  so  called  effective-modulus  theory,  the  volume  averaged 
stresses  are  related  to  volume  averaged  strains  by  means  of  effective  elastic  moduli.  Implicit 
in  the  definition  of  the  anisotropic  material  properties  is  the  existence  of  uniform 
macroscopic  state  of  stress  (Pagano,  1974).  Since  the  size  of  the  representative  element  may 
be  larger  or  comparable  to  the  region  where  the  stress  concentrations  are  dominant,  it  is  hard 
to  justify  the  validity  of  such  results.  Further  because  the  microstructure  is  ignored,  the 
detailed  localized  influence  of  the  geometric  features  simply  does  not  exist  in  these  models. 
Thus,  the  resulting  analysis  is  empirical  and  phenomenological.  Although  the  final  equations 
may  have  some  predictive  value  and  may  be  used  for  engineering  design,  they  do  not  give 
insight  into  any  microstructural  aspect  of  failure. 

From  the  existing  analytical  approaches,  it  is  clear  that  micromechanics  approach  alone 
will  not  explain  the  failure  process  of  the  laminates  simply  because  the  mutual  interaction 
between  the  micro-stresses  and  interlaminar  stresses  in  the  failure  process  has  been  totally 
neglected.  While  in  the  macromechanics  approach,  although  the  macroscopic  or  overall 
constitutive  descriptions  are  developed  from  composite  microstructure  in  terms  of  the 
volume  fraction,  the  shape,  and  the  interface  conditions  of  the  constituents,  the  constitutive 
relations  are  independent  of  the  scale  of  the  microstructure.  Further,  the  effective-modulus 
theory,  in  principle,  only  applies  to  macroscopically  uniform  fields.  Therefore,  the  stress 
fields  near  the  high  stress  gradient  regions  using  the  classical  approximation  are  questionable 
(Pagano  and  Rybicki,  1974  and  Rybicki  and  Pagano,  1975).  Since  the  details  of  the  generally 
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complex,  strongly  heterogeneous  microstructures  are  not  considered  directly,  the 
conventional  macroscopic  theory  would  inevitably  involve  erroneous  theoretical  predictions, 
and  from  which  precise  information  of  failure  in  the  micro-level  would  not  be  possible  to 
elicit. 


Apparently,  in  regions  of  macroscopically  steep  stress  (or  strain)  gradients  such  as  free 
edges  or  holes,  the  conventional  anisotropic  elasticity  theory  utilizing  the  effective-modulus 
concept  in  the  constitutive  relations  does  not  preserve  the  essential  variation  of  the 
microstress  distribution  through  the  unit-cell  in  the  micro-level.  Indeed,  rapid  change  of 
three-dimensional  stress  state  contradicts  the  underlying  assumption  of  the  macroscopically 
uniform  state  of  stress  from  which  the  effective-modulus  theory  has  been  derived.  As  a 
consequence,  at  high  stress  gradient  regions  the  conventional  approach  ceases  to  provide  true 
representation  of  physical  reality.  In  order  to  represent  the  effect  of  the  microstress  variation 
through  the  unit  cell  and  properly  capture  the  meaningful  macroscopic  steep  stress  gradient 
fields,  one  must  retain  the  volume  average  of  the  microstress  distribution  as  the  ‘stress’ 
acting  on  the  cell  but  also  the  moment  of  the  microstress  distribution,  termed  couple  stress, 
on  the  element. 

Numerical  Analysis 

In  this  work,  we  consider  the  distribution  of  the  phase  material  stresses  within  a 
simulated  free  edge  specimen.  The  fibrous  layers  are  all  oriented  in  the  longitudinal  or 
loading  direction,  while  the  remaining  layers  are  entirely  composed  of  matrix  material.  This 
analog  was  first  treated  by  Pagano  and  Rybicki  (1974).  Finite  element  analyses  shown  in  the 
figure  are  performed  under  generalized  plane  deformation.  In  the  present  case,  we  compute 
the  phase  stresses,  not  only  along  the  (artificial)  interfacial  boundaries  and  free  edge 
boundary,  but  along  the  entire  boundaries  of  the  unit  cells  surrounding  each  fiber.  It  is 
shown  that  these  stresses  lead  to  resultant  effects  in  which  both  the  forces  and  moments  are 
non-zero.  Thus,  as  shown  in  the  earlier  work,  we  demonstrate  that  the  effective  modulus 
solution,  in  which  the  composite  is  modeled  by  a  system  of  layers,  each  of  which  is  assumed 
to  have  constant  elastic  moduli,  is  not  valid.  This  may  present  significant  implications  on  the 
prediction  of  initial  damage  in  the  body  as  both  the  location  and  character  of  the  initial 
damage  differs  in  the  two  problems.  This,  in  turn,  is  shown  by  examining  the  interfacial 
stresses,  which  may  precipitate  initial  cracking  in  various  ways  depending  on  the  stacking 
sequence,  character  of  applied  loading,  and  constituent  material  properties.  The  severity  of 
the  errors  in  the  usual  effective  modulus  (ply  elasticity)  approach  are  evaluated  by 
conducting  parametric  studies  of  the  material  and  geometric  variables.  Consequences  on  the 
prediction  of  initial  failure  by  incorporating  thermal  residual  stresses  are  considered. 
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V.  Interactions/Transitions 

The  principal  investigator  has  been  working  closely  with  Dr.  N.  J.  Pagano,  Materials 
Laboratory  at  WPAFB,  on  the  technical  issues  of  microcracking  of  brittle  matrix  composites. 
The  Materials  Laboratory  has  being  tested  several  simulated  free  edge  specimens  under 
tension  normal  to  the  fiber  direction.  Failure  scenarios  of  the  failure  process  will  be 
quantified.  Currently,  Dr.  Yuan  has  been  trying  to  model  these  fracture  processes  using 
analysis  methods  and  predict  the  failure  processes  by  altering  the  material  and  interface 
strength.  The  most  dominant  singular  term  of  the  stress  field  and  associated  angular 
distributions  have  been  obtained.  Further  study  will  be  continued  to  determine  the  failure 
initiation. 
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